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^— ^ ' I begin by answering a different question, "Do we know the sphaleron rate?" and 

conclude that we do. Then I discuss a crude but purely analytic picture which 

Q.1 provides an estimate of the sphaleron rate within the context of Bodeker's effective 

fl\ • theory. The estimate, which comes surprisingly close to the numerically determined 

rf\ ' sphaleron rate, gives a physical picture of baryon number violation in the hot phase, 

and provides a conjecture of the N c dependence of the sphaleron rate in SU(JV C ) 
gauge theory. 

1 Introduction 

>. 

Originally I intended to entitle this talk, "Do we Know the Sphaleron Rate?" 

VO ■ Unfortunately the talk would then have been much too short: 

^H ■ 

G\ ! Yes. 

O 

^^ ■ Speaking seriously, it is appropriate here to remind the audience of what 

the sphaleron rate is, why we are interested in it, and how we go about calcu- 

C"| ■ lating it. As I will discuss, we know the correct effective theory for determining 

Qh the sphaleron rate; it is called Bodeker's effective theory. We also understand 

what the physics which goes into deriving that effective theory, is. Further, we 

can compute within the effective theory quite accurately. But exactly what is 

,-G ■ going on, within the effective theory, is more opaque. In the second half of the 

talk I will try to present my picture of what the physics is, which will lead to a 

purely analytic estimate for the sphaleron rate which is closer than you might 

think to the right answer. 

jrt ' Baryon number (number of baryons minus anti-baryons, essentially the 

net amount of matter) is known to be very nearly conserved, but not exactly 

conserved, within the standard model and under ordinary conditions, as shown 

long ago by t'Hoonil. The reason is the anomaly, which shows that the baryon 

number current is coupled to the SU(2) nonabelian field strength; 

where F^ v is the SU(2) weak field strength and N-p — 3 is the number of 
generations (and I have dropped an irrelevant term involving hypercharge field 

"Talk presented at "Strong and Electroweak Matter 2000", Marseille, France, June 14-17, 
2000 



strength). Further, the total derivative on the right hand side, integrated over 
spacetime in a vacuum to vacuum process, need not be zero; the topology of 
the group SU(2) shows that it can be any integer, and equals the instanton 
number of the vacuum to vacuum process. Hence baryon number is violated, 
but in vacuum to vacuum processes its violation only occurs nonperturbatively. 
Since the weak coupling is actually small, a w ~ 1/30, the instanton suppres- 
sion factor exp(— 4w/a w ) is enormous, and baryon number violation is of no 
consequence under normal conditions. 

However, the very hot conditions relevant in the very early universe did not 
constitute "normal conditions." At temperatures in excess of the electroweak 
phase transition (or crossover) temperature T ~ 80GeV, the population func- 
tions for infrared electroweak gauge bosons are large, and nonperturbative 
physics can be unsuppressecBu. This means that baryon number violation be- 
comes efficient; in fact it is ontaj polynomially, not exponentially, suppressed 
by the size of the weak couplings. 

Since the universe has a small, nonzero, and very interesting (to us!) 
abundance of baryons, understanding any cosmological process which changes 
baryon number seems well motivated, especially when it is associated with 
physics which we all believe must be there (the standard model gauge group). 
So there has been significant interest in understanding the efficiency of baryon 
number violation in the early universe. 

The sphaleron rate, the topic of this talk, is defined to be the diffusion 
constant for the quantity on the right hand side in Eq. (Fy), 

\fdt !#*>""" ; 
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r = lim lim — — = — Tm„ , (2) 

V^oo t^oo Vt 9 

with rV B the diffusion constant for baryon number. This is related to how 
fast a baryon number excess will decay, by a fluctuation dissipation relation. 
There is a free energy cost associated with having a net baryon number; 

*--H*L (3 ) 

Averaging over possible Nb with exp(— F/T) weight gives 

(Ni) = ^- , (4) 

which must be sustained by a balance between baryon number diffusion and 
baryon number decay, resulting in 

1 dN B 39 T 
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Hence, the sphaleron rate tells us how fast a baryon number abundance will 
decay. Understanding its size is the object of the rest of the talk. 

2 Bodeker's effective theory 

The context where we want to compute the sphaleron rate is the hot standard 
model above the temperature of "electroweak symmetry restoration," which 
to good accuracy meansa that we may work in Yang-Mills theory. It has been 
known for some time that the infrared behavior of Yang-Mills theory, even at 
weak coupling a w -C 1, is non-perturbativda. The easiest-way to understand 
this is to consider the thermodynamics. It is well known u that the full path 
integral for the thermodynamics of the standard model is approximated up to 
corrections suppressed by powers of a w by a 3-dimcnsional path integral, 



Z = VAi T>$ exp(-H/T) , 

H= J d 3 x i^^ + (A*) t A$ + m 2 $ t $ + A($ t $) 



(6) 



It is also well known that this 3-dimensional theory is super-renormalizable. 
That means that it has a characteristic momentum scale, in this case g 2 T. On 
scales very large compared to this, it is weakly coupled; in fact it grows more 
weakly coupled as a power of scale (until the effective theory breaks down at 
p ~ 7rT, where it is as weakly coupled as the underlying theory); conversely, 
as we look at more infrared scales it becomes more strongly coupled, until at 
the scale g 2 T its interactions are nonperturbative. Crudely speaking we could 
say that this occurs because the occupation number of a Bosonic state with 
momentum p is ~ T/p; and that the loop counting parameter is enhanced 
by one power of the occupation number. This also suggests a nice way of 
understanding the physics which occurs at strong coupling- large occupation 
number is the classical field correspondence limit, and Eq. (g|) is the partition 
function of the classical field theory at finite temperature. n 

This correspondence led Grigoriev and Rubakov to conjecture 13 that the 
right effective theory to determine the sphaleron rate is the classical field the- 
ory, regulated for instance on a lattice. We see that this cannot be true in 
Fig. [j]; the dependence on the lattice spacing is nontrivial and does not go 
away at suitably small spacing. This physics was predicted, before the data 
for the figure were taken, by Arnold, Son, and Yaffil. They argued correctly 
that "hard" (momentum scale p ~ ttT) physics is important to the dynamics 
of the infrared gauge fields, which are responsible for Chern-Simons number 
diffusion. Though not originally presented this way, their argument boils down 
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Figure 1: V for pure lattice classical Yang-Mills theory against lattica-ispacing. It does not 
show rapid convergence to a nonzero limit at zero spacing. Data from|3. 



to fairly well known plasma physics. A plasma-abelian or nonabelian-is highly 
conducting. In a conducting medium, magnetic fields evolve slowly, because 
electric fields interact strongly with the plasma. In particular, the time evolu- 
tion of a magnetic field is related to the magnitude of an electric field by (the 
nonabclian) Faraday's law. But an electric field is quickly shorted out, because 
it generates a current. The nonabelian Ampere's law reads 

E + DxB=j. (7) 

The greater the strength of screening in the medium, the greater the im- 
portance of the current term, which because it is a conducting medium is of 
form ~ j Edt, and causes the E term to "undo" past electric fields. More 
accurately, the current is 

j a (x, t) = ml J ^ v Jdy f/ ab (x, t; x-vy, t-y) v-E 6 (x-vy, t-y) , (8) 

with v a unit vector and its integral dtt v , taken over the unit sphere. Here 
U ab is the straight line, adjoint parallel transporter with initial and final 
points shown, which renders the equation gauge covariant. Together, Eq. (R) 
and Eq... (M) give classical Yang-Mills theory with added hard thermal loop 
effectaij'tjr (One must also add a noise term, which has a form completely 
determined by equilibrium thermodynamics; the noise correlator is nonlocal 



and we will not present it, se&3't3.) Arnold, Son, and YafTe advocated consid- 
ering this effective theory, and argued that the sphaleron rate should scale as 
T oc 1/mf). 

Two numerical implementations of this effective theory existoli3. Both are 
quite complicated and I will not describe them here, though I will momentarily 
discuss their results for T. 

Bodeker has argued that it is possible to integrate out the scale gT from 
the effective theory just presented to arrive at a new, "more effective'^heorJi3. 
The physical content of his calculation (which is fairly complicated Ej'EJ'HLj) 
is that scatterings between the excitations carrying the current j, exchange 
the nonabelian charges of the current carriers. The mean rate for nonabelian 
charge randomization turns out to be 
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This is parametrically larger than the scale where nonperturbative physics 
happens, ~ g 2 T, if we are willing to expand in \og(mu/g 2 T). Doing so, we 
conclude that a current carrier scatters frequently compared to the time it 
takes to traverse the length scale we are interested in, ~ l/g 2 T. Therefore, to 
leading order the current can only depend on the local value of fields, and we 
may replace Eq. M) with a local expression, which must be of form 



j = aE. 



(10) 



For obvious reasons the constant a has been christened the nonabelian con- 
ductivity [or "color conductivity," if we think about the SU(-/V C ) theory rather 
than SU(2)]. It turns out one may determine a to next to leading log order; 
at this order it is 
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3.041 



(11) 



The value for the constant 3.041 is from El Further, on the scales of interest 
the term E in Eq. (|7|) can be dropped relative to the current term, yielding 
Bodeker 's effective theory, 

j = aE = DxB + e, (&(x,t)$(yJ)) = 2oT6 ab 5 ij 6(x-y)5(t-1?), (12) 

with the noise Gaussian and white, with the autocorrelator shown. 

To determine T within this effective theory requires dealing with the non- 
linear form of DxB, necessarily including scales where the appearance of the 
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Figure 2: Sphalpwpei rate in Bodeker's effective theory, two iattice implementations ofr-HTL 
effective theory EjI13i an d pure lattice theory interpreted as HTL effective theory (seer i). 



gauge field in D is as important as the derivative term. It is a hard nonlin- 
ear problem and the only controllable techniques known are lattice techniques. 
However as a lattice problem it is particularly amenable to solution. In par- 
ticular, the effectiv£_.theory has a good, well defined limit as the UV regulator 
is taken to infinity^, and it is possible to compute analytically a matching 
between a lattice implementation and the continuum theory so the first 
are quadratic in lattice spt^ingp, and to define J F^F^ topologically 
One finds numerically thaO 
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k' = [ 10.0 ±0.2 
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This is also in good agreement with the results of the HTL effective theory, as 
summarized in Fig. @. 



3 Model for the nonperturbative behavior 

I would argue that we understand pretty well why Bodeker's effective theory 
is the correct effective description of the the infrared physics behind baryon 



number violation. Wc understand why the IR behavior is classical, and why the 
length scale l/g 2 T naturally sets the scale for topology change. We understand 
why IR gauge field evolution is over-damped, in terms of plasma physics. And 
we understand why, in a nonabelian theory, a nonabelian Ohm's law becomes 
valid on the scale l/g 2 Tlog(l/g). But the number in front of Eq. (^3|) was 
computed via lattice Monte-Carlo means, which is sort of a "black box." It 
would be nice to have some physical picture of what is going on within the 
effective theory. Here I provide a crude picture, which leads to a crude analytic 
estimate of the sphaleron rate in Bodeker's effective theory. 

The point by point distribution of F^F^ is as ultraviolet sick as that of 
F 2 . So why is its diffusion constant finite and infrared dominated? The answer 
is that F^F^ has topological meaning. It is a total derivative and can be 
written as the spacetime gradient of a current, (g 2 /i2ir 2 )F^ 1 ' F^ u = d^K^. 
The associated charge is the Chern-Simons number, 

2 /• 

d 3 xK Q = A cs = ^Jd 3 x e ljk (F*Al - |/ q6c A^^) . (14) 

In vacuum the integral of the A 3 term is always an integer; so there must be 
nonvanishing fields for Acs to take on a non-integral value. And in particular, 
since FijAk has one fewer derivative than F 2 , the more concentrated the region 
where Acs is stored, the greater the energetic cost. Physically, this means that 
a fluctuation which creates a small amount of Acs leaves an imprint or "strain" 
behind in the fields, so there is gradient energy present until either the fields 
relax in a way which undoes the Acs generated, or enough more is added to 
get up to the next integer. This is what prevents there from being any UV 
divergences-in fact, any unsuppressed UV contribution at all-in T. 

The other piece of physics we know about the symmetric electroweak 
phase, or pure 3-D Yang-Mills theory, is that there is magnetic screening. 
Colored degrees of freedom interact with each other exclusively through gra- 
dient interactions, so a colored field at x feels the fields at y only through a 
parallel transporter from y to x. In fact, the field equations mean that the 
fields at x will generally evolve under the influence of y through some average 
over paths for parallel transportation. 

For x and y of small separation <C 1/ g 2 T, two typical paths between them 
will have approximately the same parallel transporter. Therefore, the interac- 
tions via the different paths will add coherently. But for x and y separated 
by <; 1/ ' g 2 T ', a generic pair of paths give completely different parallel trans- 
portations. (Equivalently, the mean trace of large Wilson loops goes to zero.) 
Roughly, there is a length scale / m ag, called the magnetic screening length, 
beyond which parallel transportation differs by 0(1) on change of path. This 



means that the influence of fields at y on fields at x add incoherently over paths 
and will approximately cancel. A field cannot "see" colored objects more than 
l mag away. This arises because of large nonabelian fields in between the two 
points, which are responsible for making the parallel transportations differ. 

The physical picture I argue for is that magnetic screening renders the 
"gradient energy" argument for keeping Ncs from changing, ineffective for 
field fluctuations on scales <; Z mag . Hence, Ncs will diffuse due to F^F^ 
from fields at the scale I. These fields will evolve diffusively; but fields on less 
IR scales will not contribute to the diffusion of iVcs • 

Within this picture, we can make a crude calculation of T, as follows. 
Consider some large volume V as & collection of boxes l mag on a side. Within 
each box, the most IR gauge field degrees of freedom diffuse, and all higher 
wave number modes do not contribute; and the contribution of different boxes 
are independent. Writing F^F^ = AEfBf, Eq. (§) becomes 

^/^"■'''^'■■'''tV (15) 




8tt 2 

where the subscript on the magnetic field means that only the most infrared 
component is kept. The sum over directions gives the two independent trans- 
verse polarizations, giving 

^L_ Y, Jdtidt 2 dx 1 dx 2 (E a (yL l ,ti)E b (x 2 ,t 2 )B? R (x l ,t l )B b K (x 2 ,t 2 )). 

boxes 

(16) 
We determine the E field correlator from Eq. ( |f2| ) . Our argument about scales 
says that, if we consider only the most infrared component of the magnetic field, 
then we can neglect the gradient term, D x B, in evaluating the E correlators. 
In this case E and B are uncorrelated, and the E correlators come only from 
the noise, 

IT 

(E a ( Xl ,h) E b {x 2 ,t 2 )) = —5 ab 5{x l - x 2 )8{t l - t 2 ) . (17) 

a 

(For higher frequency excitations the D x B term is what will ensure no con- 
tribution.) This performs one space and one time integral: 

fdxi dx 2 dh dt 2 ( E a (x u h) E\x 2 ,t 2 ) B? K (xi,h) B b lK {x 2l t 2 ) ) 

IT f 

= —5 ab dx 1 dx 2 dt 1 dt 2 5(x 1 -x 2 )5(t 1 -t 2 ) (B^ R (x u t 1 )B^ R (x 2 ,t 2 )) 

= ^S ab f dhjdx {B^{x,h)B h IK {x,t 1 )) . (18) 



This integral can be estimated by equipartition, 



J dh [d 3 x (B? R (x,ti)B% R (x,ti)) = TS ab J dh = Ttt 
Jo J Jo 



-a b 



(19) 



Probably there is more power in the infrared magnetic fields than equipartition 
would estimate; but this is only one point at which our analysis is crude. 

The group indices have reduced to S ab S ab = N%—1. The volume of a box 
is ?magi so m a total volume V there are V/l^ lag boxes. Putting everything 
together, we get 



r = 



4(7V c 2 -l)g 4 r 2 
which, on substituting Eq. (O) and re-arranging, gives 
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3.041 



(21) 



It remains to estimate Z mag , which is a thermodynamic question and can 
be answered within the 3-D Euclidean theory, Eq. (g) . The only natural length 
scale of hot Yang-Mills theory at weak coupling is l mag ~ l/g 2 T, and at large 
N c the N c dependence must be l ma , g ~ l/N c g 2 T. The only question is what the 
coefficient is. I estimate it by trying to find on what length scale the perturba- 
tive expansion for 3-D Yang-Mills theory (describing the thermodynamics of 
the thermal theory) has completely broken down. 

One estimate is that Z mag ~ 1/p for the momentum p with |ITr(p)| = p 2 , 
that is, where the one loop transverse self-energy computed with tree prop- 
agators is as large as the tree inverse propagator. In Landau gauge (strict 
Coulomb gauge from the 3+1 dimensional perspective) this is 



II(p, one loop, tree propagators) 



_ HN c g 2 T 

64|p| 



(dijP 2 - PiPj) , (22) 



so l mag ~ 64:/(llN c g 2 T). This estimate is gauge dependent; in Feynman gauge 
the 11 becomes 14. However, since unequal separation correlators are maximal 
in Landau gauge, it seems the appropriate gauge for estimating where those 
correlators break down. 



of the "magnetic mass" 
Without endorsing the 



Alternately, we can take / mag as the reciproc 
found self-consistently by solving a gap equation E^ 
efficacy of gap resummed perturbation theory to make useful predictions, one 
can nevertheless say that the mass found gives a characteristic scale where 
higher order terms in the perturbative expansion are of order the leading ones. 
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Figure 3: Volume dependence of T (pure lattice theory, g 2 aT = 1/2) illustrates that V turns 



A two loop gap equation calculation E3 gives Z mag = l/rn mag = l/0.17N c g 2 T 
which agrees with the self-energy estimate!] 

One check on this estimate of Z mag is to look at lattice data in finite volume, 
and see how large a lattice has to be before T "turns on." We expect that it 
should have to be about 2Z mag across, because we use periodic boundary con- 
ditions on the lattice. Fig. indicates that our estimate of Z mag is reasonable. 
In particular, the inverse glueball mass < 1/ ' g 2 T is not the appropriate scale. 

Using the estimate for Z mag , we get 



T- .247V c 3 (iV c 2 -l) 



N c g 2 T 2 



•d 



aiT 4 



l 0g 1^ + 3.041 
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(23) 



Plugging in iV c = 2 gives V = 11.5(g 2 T 2 /m 2 D )al ! T 4 \og{. . .). The 11.5 should 
be compared with the actual number, 10.0, presented earlier. The quality of 
the agreement is probably fortuitous. 

4 Conclusion 

After reviewing the current status of what the sphaleron rate is, and what 
the right effective theory for its determination is (Bodeker's effective theory), 

6 The computation^ is in SU(2), but since all graphs at one and two loops are planar, the 
N c dependence must be as shown. 
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I have proposed a physical picture of what is going on in the strongly coupled 
IR physics which determines the rate. My picture is that, on scales longer 
than Z mag , magnetic screening prevents "any record of being kept" of what 
past changes have occurred to Ncs', therefore, very IR gauge fields evolve 
diffusively. I have presented a crude estimate of this diffusive evolution and of 
( mag , and find that the resulting sphaleron rate T is surprisingly close to the 
value determined in lattice simulations. 

The model here should not be taken as an accurate quantitative predictor 
of T; rather it should be viewed as a parametric estimate where an effort 
has been made to include factors of 2 or 7r where possible. It also gives an 
interesting physical picture of Nq§ diffusion in Yang-Mills theory; rather than 
occasional, integer changes, Nq$ diffuses from the incoherent accumulation of 
many contributions individually much smaller than an integer. 

A more robust prediction of this analysis is for the dependence of T on 
N c . Of course the scaling behavior I ~ l/N c g 2 T has 0(1/ N 2 ) corrections, so 
the N c dependence in Eq. (p3J) cannot be relied on absolutely at small N c . 
However, it would be interesting to numerically investigate T within Bodeker's 
effective theory for a few larger N c and see how welLthe prediction holds. In 
particular the case of SU(3) is physically interesting^! 

In conclusion: Do we understand the sphaleron rate? I think, perhaps, 
that we do. 
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